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ABSTRACT 


In this paper we first derive two normal form constructions 
for cellular automata to transform any given (one-dimensional) cel- 
lular automaton into one which is one-way and/or totalistic. An 
encoding of this restricted type of automaton together with any ini- 
tial configuration becomes the input of our small universal cellular 
automaton, using only 14 states. This improves well-known results 
obtained by simulation of small universal Turing machines and also 
some recent results on universal totalistic cellular automata. 


1. Introduction 


Interest in cellular automata (CA) has been renewed since their application 
to the study of complex systems [19-21]. In this context the universality of CA 
was discussed in [21] and open problems about universal CA with a small number 


of states were stated in [22]. We give some results here. 


A cellular automaton (CA) is said to be universal if it can simulate every 
Turing Machine or, even stronger, if it can simulate every CA of the same dimen- 
sion. The first universal cellular automaton was given in the famous work of J. 
von Neumann on the simulation of self-reproduction. He gave a 2-dimensional 
universal and self-reproducing CA with 29 states. This was improved to 20 states 
in [1]. 

Also well known is the work on small universal Turing Machines (cf. [11, 


12, 15]). In this case the goal is to minimize simultaneously the number of states 


and the number of tape symbols. In [12] it is shown that there is a universal 
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Turing Machine with either 4 symbols and 7 states or with 6 symbols and 6 
states. Smith [16] has shown that any Turing Machine with m symbols and n 
states can be simulated by a one-dimensional CA whose uniform cell has m 4-2n 
states. Thus either one of the universal TM from [12] yields a universal CA with 
4+2.7 = 6-+2.6 = 18 states. Here we improve this result and show that there is а 
universal one-dimensional CA with 14 states. We believe that some more effort 
in the detailed “low-level programming” of our basic strategy can decrease this 
by 2 or 3 states. However, the minimal number of states of а universal one- 
dimensional CA could still be much smaller, since already CA with 3 or 4 states 


show astonishingly complex behaviour. 


The well known "Game of Life" is a "semi-totalistic" CA which means that 
the next state of a cell only depends on its own current state and the sum of the 
states of its neighbors. S. Wolfram also introduced an even more restricted type 
of CA, called totalistic. The next state of a cell of such CA depends only on the 
sum of all the states in its neighborhood, including its own. D. Gordon [9] has 
shown that totalistic one-dimensional CA can simulate every Turing Machine and 
we will strengthen this and show that every one-dimensional CA can be simulated 
by а totalistic one. We will use this result in the construction of our small 
universal CA. We will also use the result that every one-dimensional CA can be 
simulated by a one-way (unidirectional) one-dimensional CA [6,8,18] and that 


one-way CA are equivalent to trellis automata [4]. 


2. Preliminaries 


As we will consider only one-dimensional and homogeneous cellular auto- 
mata in this paper we will restrict ourselves to this special case in the following 


definitions. The more general terminology can be found in [5] or [17]. 


Intuitively, a cellular automaton consists of a doubly infinite array of cells. 
All cells are identical copies of one single finite automaton. The local transition 
function of each cell only depends on the actual states of its left and right neigh- 
bor and itself. Thus, à computation of a cellular automaton can be defined in the 
straightforward way as the synchronous application of the local transition func- 


tion at each cell. The set of states always contains a so-called quiescent state q 
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with the property: if a cell and its left and right neighbors are quiescent at time t 
then this cell is quiescent at time t+1. As we assume finite input for the cellular 
automaton this implies that there is a finite number of nonquiescent cells in the 


initial configuration and in every subsequent configuration as well. 


Definition 2.1 A (one-dimensional, homogeneous) cellular automaton is a triple 
А = (Q, d, q), where Q is a finite set of states, d is the local transition function, 
d: QXQXQ — Q where the arguments of d are used in the following meaning 
d(state of left neighbor, own state, state of right neighbor), and q in Q is the 
quiescent state, i.e. it holds d(q, q, q) = q. 


A configuration of A is a mapping C : Z — Q, where Z denotes the set of 
integer numbers such that C(i) = q (quiescent state) for all but finitely many 15. 


The set of nonquiescent cells of a configuration C is called the support of C. 


А computation of A now consists of a sequence of configurations 
CC ,C5...,C,,..., Where Cg is the initial configuration and each configuration 
С; is generated by simultaneous invocation of the transition function d for all 


cells of A with states as given by C;. 


3. Simulation by one-way automata 


For our construction of a universal cellular automaton in subsequent 


chapters we will need an automaton which is totalistic and one-way. 


A cellular automaton А with set of states Q and transition function 
d:QxQxQ — Q is called one-way, if d only depends on the state of the own 
cell and the state of its right neighbor cell. Thus we can write the transition 
function of а one-way cellular automaton in the form e: QXQ — Q with the 
convention that the arguments of e consist of the states of the own cell and those 


of its right neighbor. 


In this section we will outline а transformation of an arbitrary cellular auto- 
maton to a one-way cellular automaton in order to specify bounds for the 


increase in the number of states and in time-steps needed in the simulation. 


In [18] а similar technique was used for the case of real- time cellular auto- 


mata, more general cases were considered in [6] and [8]. 
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The construction of the one-way automaton can be described briefly by: 
1. merging of each state with the state of its right neighbor 
2. shifting to the left during the state transition. 


Let 8,,55,..,8, be states in a configuration, such that s, Æq and let 


to, £1,...,¢44, be the states after one transition as shown below 


q q q 8, 82 83 `` 8k q q q 

q q to ti t2 ts coc t aq 9 
Then it takes 2 time-steps in our simulating one-way automaton to produce the 
transition | 

q q q Sí 829 83 ttt 8k q q q 

q to ti t2 ts Uo f. Q q q 


For the given cellular automaton A let Q be the set of states and 


4: QXQXQ — Q the transition function. 
As defined above, in the simulating one-way cellular automaton Æ the tran- 
sition function d' will consider only the states of its own cell and its right neigh- 


bor cell. 
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Let Q! = QU QXQ and define d' : «О! + Q' as follows 
d(u,v) = uv for (u,v) я (4,4) 
d(q,q) = q 
(ио, vw) = d(u,v,w) 
d'(q,qu) = d(q,a,u) 
d'(wq,q) = d(w,q,q) 


for all u,v,w in Q. 


Thus, the transition shown above is completed in two time-steps by A in 


the following way 


q q 81 82 83 Sk q q q 
q 981 8182 8283 Sk q q q q 
to t, to cf. fu Q q q 


The details of this construction are straightforward, so we can state the fol- 


lowing theorem: 


Theorem 3.1 For every cellular automaton А with k states there exists a one- 
way cellular automaton А which simulates A twice slower and A needs at most 


k? + k states. 


In the above construction of the simulating one-way cellular automaton one 
can delay the given transitions by "aging" the states іп Q (eg. а, a’, a", а",...) to 


obtain slower expansion of the nonquiescent cells. 
Therefore it is clear that the following corollary holds: 


Corollary 3.2 For every cellular automaton А with k states and for every 
m>1 there exists a one-way cellular automaton A which simulates А(т-+1)- 


times slower and A needs at most К? + mk states. 


Later, we will have to apply this corollary for m=3 in the construction of a 
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universal one-way cellular automaton. 


4. Simulation by totalistic automata 


In this section we will show that each cellular automaton can be simulated 
(without loss of time) by a cellular automaton which has a totalistic transition 
function and uses up to four times as many states as the original automaton. 
Definition 4.1 A cellular automaton A with set of states Q and transition func- 
tion f : QXQXQ — Q is called totalistic, И Q C М (non-negative) and there 
exists a function f! : N — М such that f(z,y,z) = f'(z--y--z) for all z,y,z in Q. 

The transformation of an arbitrary cellular automaton whose states are 
identified with non-negative integers is now accomplished by a cyclic "coloring" of 
cells. We use four different powers of a basis such that in the number representa- 
tion of the sum of three neighboring states left neighbor, right neighbor and own 


cell are still identifiable by the position of a missing entry. 


Consider e.g. a cellular automaton A with set of states Q = (1,2,...,n) and 
transition function 4: ОХОХО — Q and let the figure below depict part of а 


configuration of А. 


а и ме ESL ERAS ESL SE 


Let B = п+1 be the basis for the colouring factors 10, 100 and 1000 (in B-ary 


notation). Then the above configuration changes to the following: 


—> — — — —» 
1008 10e., 10008, 1083 
——— < —— <—— 


Thus our new set of states 1s 


Q! = (sm |s ЕО, т Є (1,10,100,1000)) . 
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So @ contains 4n states and we can define now the (partial) totalistic transition 


function d! : N — Q! as follows 
d'(xyz) = 10d(z,y,z) 
d'(xyz0) = 100d(z,y,z) 
d'(yzx) = 1000d(z,y,z) 


d(z0zy) = d(z,y,2) 


for all z,y,z in Q. 


Again zyz, zyzO, yzOz, zO0xy, 10, 100, 1000 are to be interpreted as 


numbers in the B-ary system. 


It is now straightforward, that the cellular automaton A with set of states 
Q' and totalistic transition function d' correctly simulates A without loss of time 


and it therefore holds 


Theorem 4.2 For every cellular automaton A there exists a totalistic cellular 
automaton А which simulates А without loss of time and has at most four times 


as many states as A. 


If only one-way cellular automata are considered, the above construction 
can be simplified using only three different powers of B for the colouring of the 


cells, whereby the one-way property of the given automaton is preserved. 


Corollary 4.3 For every one-way cellular automaton A there exists a one-way 
totalistic cellular automaton A which simulates А without loss of time and has at 


most three times as many states as А. 


For the more general case of an arbitrary (regular) systolic network this 
technique of coloring its underlying graph is used in [23] to show that actually 


every systolic network can be transformed into a totalistic one. 
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5. Informal description of construction 


In the previous section we have shown that any cellular automaton can be 


simulated by one which is one-way and totalistic and such that the expansion of 


the nonquiescent part to the left can be delayed by any constant factor. We 


choose m —3 to achieve expansion to the left at most at half speed. 


Our universal cellular automaton U will simulate any given one-way totalis- 


tic cellular automaton А ( which expands at most at half speed to the left) with 


any given initial configuration. We will encode this pair as an initial configuration 


of the universal cellular automaton U. A constant number (depending only on 


the number of states of А) of steps of U is needed to simulate one step of A. 


The time-space diagram (unrolling) of the given one-way automaton then looks 


like 
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This unrolling is obviously equivalent to the following 


О 


Ni UP GUN 


ААА ТАТА ТА $4. 
x» 
\ 
ИМУУУМУ N 


PLS LL oT 
OY 


which will be the basis of our construction of the universal cellular automaton U. 
Note that the computations which took place in one cell of the first automaton. 
therefore being depicted in vertical lines in the first figure, are now shifting to 


the right in each step. 


Example: Let A be a totalistic one-way cellular automaton with set of states 


Q = {0,1,2,3,4} and transition function d given by 


i JO 1 2 3 4 5 6 
di)|o o 1 3 2 3 8 


Since A is one-way the nonquiescent states expand to the right only. An exam- 


ple of а computation of А follows. 


VAVAVAVAVANAN 
ATATAVATAYAY 
ЛАЛАЛАЛ 
`\АЛАЛАЛА 
ГАУАТАТАТАТАТАТА 
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In the simulation of the automaton А each of the nonquiescent cells of A will be 
encoded as a block of constant size of cells in U. This block contains the given 
transition-table for the cells of A, and by the marking of the appropriate table- 
entry also the encoding of the current state of a cell, ie. if a transition 
d(i4-j) = k takes place in a cell of A, in the corresponding block the (+7) 


table-entry - which holds an encoding of the number k-becomes marked. 


So in our automaton U the main actions to simulate one time-step of the given 


automaton A will be grouped into three phases 


phase 1: send information about current state to left neighbor block 
phase 2: send information about current state to right neighbor block 
phase 3: process information to look up next state. 


6. Outline of construction 


As shown above the next states are always computed between the two cells 
whose states are to be considered for the table-lookup. Therefore it is convenient 
to use alternatively “activated” blocks and “empty” blocks and to change the 
roles of both in the odd and even steps (beats) of the computation. This is 
achieved by sending messenger signals over these blocks where the signals have 


to complete various tasks during the three phases mentioned above. 


phase 1: start from the right end of an active block and find the table-entry 
representing the current state; activate a number of signals according 
to the current state which are moving to the empty block to the left; 


continue to move to the left end of the block. 


phase 2: bounce at the left end of the block and find once more the table- 
entry representing the current state; accordingly activate signals 
which are now moving to the empty block to the right; while clearing 
all temporary marks continue to move to the right end of the block, 


leaving an empty block behind. 
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phase 3: cross over to the block to the right; while moving to the right end of 
this block process the signals which have been accumulated here dur- 


ing the previous phases; bounce at the right end. 


The flowing of messenger signals over some adjacent blocks during the phases of 
[ š : 
i two consecutive beats can be depicted as follows: 


block š block i41 block i42 block $43 block $44 


Furthermore we have to take care that after each beat a new block at the right- 
hand end of the nonquiescent part can be activated to simulate a possible expan- 
sion to the right. This is achieved by copying (at least) one complete block to the 
right during each beat and by activating one new block which starts in the 


; encoding of the quiescent state. 


By our construction we can assume that there is no expansion to the left and so 
we only need to activate the leftmost block at each second beat. This is realized 


by а special block and messenger signals. 
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Details of these constructions and encodings, the transition-table for U and 


examples are found in the appendix. 


Theorem 6.1: There exists a universal cellular automaton U with 14 states 
which can simulate any given totalistic one-way cellular automaton A with any 


initial configuration. 


In [16] it is shown (Theorem 4) that any Turing machine T with m tape-symbols 
and n states can be simulated by a cellular automaton with m+2n states. Thus 
either of the universal Turing machines with 6 symbols and 6 states or with 4 
symbols and 7 states (cf. [12]) yields a universal cellular automaton with 18 


states. 


We can now apply the normal form constructions of our previous sections 
to the universal cellular automaton U of Theorem 6.1 to generate universal cellu- 


lar automata which are also one-way and/ or totalistic. 


Corollary 6.2: For every universal cellular automaton U with n states there 


exists a universal one-way cellular automaton U, with n? + n states. 


Corollary 6.3: For every universal cellular automaton U with n states there 


exists a universal totalistic cellular automaton U, with 4n states. 


Corollary 6.4: For every universal cellular automaton U with n states there 


exists a universal one-way totalistic cellular automaton U, with 3(n? + n) states. 


7. Appendix 


For the encoding of the given transition-table as a block of cells in U we 
can assume w.l.o.g. that for the given totalistic one-way cellular automaton A the 
set of states Q is defined as Q = (3,4,5,...,k), where 3 is identified with the quies- 


cent state in A. 


Each block consists now of a sequence of 2k+1 sub-blocks of the form 


#2z000...000111...111 
kama comand and a block-endmarker 4, where z is in (0, 1}, 
i—limes j—times 


¿+j = 2k--1 and if in the given transition function d we have 4(т)=п then in 
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the mth sub-block of each block it holds j7=n. For all other sub-blocks we have 
2j—3. Thus a block consists of (2k+2)° cells. 


If in a cell of the original one-way totalistic automaton a transition d(r)=s 
takes place then in the corresponding block in U exactly in the first r sub-blocks 
from the left this z equals 1 and in all other sub-blocks z=0. 

The computation of the next transition in the simulating automaton U then 
takes 3(2k--2)? + 7(2k--1) time steps which constitute one "beat" of U. 


Example: Let О = {3,4,5} be the set of states and the transition-function d be 


defined as 
i 6 7 S8 9 10 
di)|3 4 4 5 4 


Then a transition by 4(9)=5 corresponds to the following situation in the block 
at the beginning of the beat: 


#1000000001 11 #1000000001 1 13 100000000: 1 1 #1000000001 1 1 #1000000001 1 1 # 10000000011 1 
#100000001 1 113€ 100000001 1 1 1 #10000001 1 1 1 1400000000 1 1 1 1 #0000000001 aad 


Here aa at the end of the block shows that this block is activated and that the 


new state 5 has to be sent to left and right neigbors. 
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It can also be seen easily that the transition-function has been extended to 


i |1 2 3 4 5 6 7 8 9 10 H 
di)|3 3 3 3 3 3 4 4 5 4 3 


An inactive, "empty" block then looks like 


#0000000001 1 1 44000000000 1 1 1 #0000000001 1 1 440000000001 1 1 #0000000001: 1 1 00000000001 1 1 
#000000001 1 1 14000000001 1 1 1 #00000001 1 1 1 184000000001 11 1 #0000000001 1 1 # 


If a block is activated at the beginning of a beat, his right neighbor block will 


become activated at the next beat. 


In order to be able to expand, it is therefore necessary to start a new 
activation from the left end at each second beat. This is done simply by sending 
a synchronizing signal up and down over a sufficiently long sequence of 1’s. The 
first block after this synchronizing left part contains only table-entries represent- 
ing the quiescent state. So each second beat the "quiescent block" at the left end 


is activated and no expansion to the left can occur. 


At the right end copying of a block is performed by pushing an encoded 
form of a block to the right and simultaneously leaving an empty block behind, 
which will be activated later as representing the original quiescent state. During 
one beat a little more than a complete block is created at the right, but since . 
activation of blocks is initiated from the left - synchronized with the beat - this 
does nof affect the simulation of the nonquiescent part of the original totalistic 


one-way automaton. 


Example: The following three snapshots show the simulation of two transitions 


\/\ 
/NPN 
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where the transition-table is given as ` 


In this case one beat of U consists of 363 single transition-steps. The following 
are the initial configuration of U, the configuration after the first beat (363 


steps), and after the second beat (726 steps). 


q399111111111111111111111111111211111111111111111111112121111111111111111111111112 212121113 HE 2L 2LEI EL LII 
1111111111111111111111121112111111111111111111111111111112111111111111111111111] 11 EL ELEEL E 2211 2E IL 21] 
111111111111111111111111111111111112111111111111111111111111111111111111211111111111111111111111111 11 
111111111111111111111111221111111111111111111111111111111 1H LLL E LL 

1400011111111#000001 113400000001 1 1 40000000 1 1 14400000001 113400000001 1 14400000001 1 14&00000001114$«00000001 ия 
# 10000001 1 1 #10000001 1 1 # 10000001 1 1 #10000001 1 1 #10000001 1 1# 10000001 1 1440000000 1 1 1 #00000001 1 1 #00000001 aa sf 
#0000000 1 1 1 #00000001 1 1 #00000001 1 1 #00000001 1 1 #0000001 1 1 1 #00000001 11#0000001 111 #0000001 1 1 1 #00000001 1 1# 
# 1000000 1 1 I #10000001 1 1 #10000001 1 1 #10000001 1 1 # 1000001 1 11440000000 1 1 1 #0000001 1 11440000001 111440000000 i aa# 
$$ 10000001 1 1# 1000000 1 1 1 # 10000001 1 1 #10000001 11 #1000001 1 1 1 #10000001: 1 1 #1000001 1 1 1 #10000011 1 1 #10000001 1 1# 
#10000001 1 1# 1000000/g 1g0h0h0h0g 1g 1hOhhhhhggg | hOhOhOg 1g 1gOhhOhOhggaghhhOhOg ig 1g0hOhOh 1gg 1qqqqqqaqaqqqq 
qqaqqqqqqqqqqqqqqqqqqqqqqqqqqqqqqqqqqqqqqqqqqqqqqqqqqqqqqqqqqqqqqqqqqqqqqqqqqqqqqqqqqqqqqqqqqqqqqqqqq 


qqqq111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111 
12111111111111111211111111111111111111111111111111111111311311111111111111111111111 2211111111111 11111 
11111111111111111111111111111111111111111111g1111111111111111111111111111111111111111111111111111111 
11111111111111111111111111111111111111111111111111111111111111111944 

1440001111111144000001113$400000001 11#1000000111#10000001 1 1 4 10000001 1 1 #10000001 1 1 # 10000001 1 1 #100000011 1# 
#0000000 I 1 1 #00000001 1 1 #00000001 1 1 #00000001 1 1 440000000 1 114400000001 1 1 #00000001 1 194000000011 1 #000000011 1# 
# 10000001 11 #10000001 1 1 4 10000001 1 144 10000001 1 #1000001 1 1 1 #10000001 ! 1 #1000001 1 1 1 #0000001 1 1 1 #0000000 1aa# 
#00000001 1 1 #00000001 1 1 #00000001 1 14400000001 1 1 #0000001 1 114400000001 1 i #0000001 1 1 1440000001 1 1 1 #000000011 1# 
#10000001 1 1 #10000001 1 1 #10000001 1 1 #10000001 1 1 #0000001 1 1 1 #00000001 1 1 #0000001 1 1 1 #0000001 11 1 #00000001aa# 
#10000001 1 1 #10000001 1 1 #10000001 1 1 # 10000001 1 1# 1000001 1 1 1 #10000001 1 1# 1000001 1 1 1 #1000001 1 1 1 #100000011 1# 
#10000001 1 1 ¥ 10000001 1 1 # 10000001 1 1 #bjchOhbg lgghhhhhh1gig1hOhOhOh 1g ighhhbhhg1g1gOhOhCh 1g1gThhhhhhhg]1g] 
В110)99999999999999999999999999999999999999999999999999999999999999999999999999999999999999999999999 


інши 
11111111111111111111212111112222E222211222220100 20211122202 22222EELLLLLEELLL ELE LLL LLL LLLI 
J111111111111111111111111111111111111111111111111111311111111111111111111111111111111111111111111311 
пинин 
1400011111111#00000111#00000001 1 1 #00000001 1 1 440000000 1 1 1 #00000001 1 14800000001 1 1 #00000001 1 1400000001 1 14k 
# 10000001 1 1 #10000001 114410000001 1 1 #10000001 1 1# 10000001 1 1 ¥ 10000001 1 1440000000 1 1 1 #00000001 1 1 #00000001aa# 
#00000001 1 1 #00000001 1 1 #00000001 1 1 #00000001 1 1 440000001 1 1 1 #00000001 1 1 #0000001 1 1 140000001 1 1 1 #000000011 1 ¥ 
# 1000000 1 1 144 10000001 1 1 # 10000001 1 1# 10000001 1 1 #1000001 1 1 1 #10000001 1 1 #1000001 1 1 1 #00000011 1 1400000001 aa# 
#00000001 1 100000001 1 1 440000000 1 1 1 #00000001 1160000001 1 1 1 #00000001 1 1 #0000001 1 1 1440000001 1 1 1 #00000001 1 1 # 
W 10000001 1 1 #10000001 1 1 #10000001 1 1 #00000001 1 1 #0000001 1 1 1 #00000001 1 1 #0000001 1 1 1 #0000001 11 1 440000000184 
#10000001 1 1 #10000001 1 1 #10000001 1 1 #10000001 1 1 #1000001 1 1 1% 10000001 1 1% 1000001 1 2 1# 1000001 1 1 1410000001 1 1# 
# 1000000 1 1 1410000001 1 1 #10000001 1 1 #10000001 1 1 #1000001 1qgghhhhohOg 1g 1 hOhOhOggggghhhOh0g1g1g0hOhohOgege 


hghhOhOh1g1gOhOhOhOgg 1qqqqqqqqqqqqqqqqqqqqqqqqqqqqqqqqqqqqqqqqqqqqqqqqqqqqqqqqqqqqqqqqqqqqqqqqqqqqqq 
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The following table of state transitions for our universal cellular automaton 


table for U the upper row shows 


lon- 


* 


In the transit 


U will need some explanations. 


the state of the considered cell and the first two columns display the states of the 


left and right neighbor cells. Blank table-entries mean that these combinations 


never occur during a simulation, thus any of the states could be entered instead 


of the blank. 


U 1 a bcd e f q h i 


= 55 
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